Sign-reversal of the in-plane resistivity anisotropy in hole-doped iron pnictides 
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The in-plane anisotropy of the electrical resistivity across the coupled orthorhombic and magnetic 
transitions of the iron pnictides has been extensively studied in the parent and electron-doped com- 
pounds. All these studies universally show that the resistivity pa across the long orthorhombic axis 
ao - along which the spins couple antiferromagnetically below the magnetic transition temperature 
- is smaller than the resistivity pb of the short orthorhombic axis bo, i- e. Pa < Pb- Here we report 
that in the hole-doped compounds Bai-ccKccFe2As2, as the doping level increases, the resistivity 
anisotropy initially becomes vanishingly small, and eventually changes sign for sufficiently large 
doping, i. e. pb < pa- This observation is in agreement with a recent theoretical prediction that 
considers the anisotropic scattering of electrons by spin- fluctuations in the orthorhombic/nematic 
state. 



INTRODUCTION 

The undoped and moderately-doped 74Fe2As2 (^=Ca, 
Sr, Ba, Eu) compounds undergo a tetragonal-to- 
orthorhombic transition at a temperature T^, which is 
either accompanied or followed by stripe-type magnetic 
order below T/v < [1 , with in-plane ordering vector 
(tt, 0) or (0,7r). Since superconductivity in these com- 
pounds is induced upon suppression of the structural and 
magnetic orders by doping or pressure, it is natural to as- 
sume that magnetic fluctuations play an important role 
in the pairing interaction. Therefore, it is of great impor- 
tance to understand not only the origin of the magnetism 
of these materials, but also how it affects the normal state 
properties. Interestingly, experiments in different types 
of materials suggest a close relationship between uncon- 
ventional superconductivity and spin-fluctuation scatter- 
ing [2. 

In the phase diagram of the iron pnictides, the mag- 
netic and structural transition lines always follow each 
other OS]. To describe the close relationship between 
these two orders, two alternative scenarios have been pro- 
posed. One scenario considers that the structural tran- 
sition is a secondary consequence of an emergent Ising- 
nematic phase, which appears in the paramagnetic state 
due to the interplay between magnetic fluctuations and 
the doubly-degenerate character of the striped magnetic 
phase . Another scenario considers that spontaneous 
orbital ordering triggers the structural transition and re- 
moves the intrinsic frustration of the striped magnetic 
state, which then appears at lower temperatures [9-12 . 

Since the nematic, orbital, and orthorhombic orders 
all break the tetragonal symmetry of the system, it is 



difficult to decide between the different scenarios only in 
the qualitative level. Thus, it is fundamental to exper- 
imentally investigate the macroscopic properties of this 
anisotropic state thoroughly. In this context, the resis- 
tivity anisotropy is an important benchmark, since sev- 
eral concurring mechanisms for anisotropic transport are 
present, such as Fermi surface reconstruction in the mag- 
netically ordered state [TSHIS], electronic reconstruction 
due to orbital order [TTJ [T2j [TBHTS] , and anisotropic scat- 
tering by spin fluctuations [19]. 

Experimentally, the resistivity anisotropy can be stud- 
ied in crystals which have been detwinned by the applica- 
tion of an uniaxial stress or strain [15 , 20 , see Ref.[21 for 
review. A very unique picture was observed in the parent 
and electron-doped materials [151 [20] ^(P^i-x TM x)2^^2 
(TM=Co, Ni, Cu): counter-intuitively, the resistivity in 
the antiferromagnetic state is larger along the direction of 
ferromagnetic chains (short orthorhombic ho direction) 
than along the direction of antiferremagnetic spin chains 
(long orthorhombic ao direction). The anisotropy re- 
mains finite well above the orthorhombic transition tem- 
perature of the unstressed sample, suggesting the pres- 
ence of a large nematic susceptibility in the tetragonal 
phase, in agreement with elastic modulus measurements 

In our initial study [15], we concluded that the 
anisotropy in the magnetically ordered state can be 
well understood using band structure calculations. De- 
spite the good agreement of this conclusion with recent 
Shubnikov-de Haas oscillations measurements [23] and 
detailed transport analysis [24 in detwinned crystals of 
undoped BaFe2As2, the authors of Ref.|22 argued that 
with this mechanism it would be very difficult to explain 
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the large in-plane anisotropy ph/ pa ^ 2 observed in Co- 
doped samples. Alternative explanations focused on the 
unequal occupation of dxz and dyz orbitals in the antifer- 
romagnetic phase as the origin of the anisotropic trans- 
port ^13, .24,. 

While these discussions focused on the effects of the re- 
constructed band structure on the resistivity anisotropy, 
the temperature dependence of Ap = — p^, and also 
optical conductivity data [25 , reveal the importance of 
anisotropic scattering as well. The resistivity anisotropy 
in all parent compounds [26 is maximum at Tg, where 
both the orthorhombic lattice distortion and the magni- 
tude of the magnetic moment are minimum. This mo- 
tivated a model where scattering by spin fluctuations, 
which are anisotropic in the nematic/orthorhombic state, 
plays a fundamental role in determining the properties of 
the resistivity anisotropy [19j. An interesting prediction 
of this model is that the resistivity anisotropy in hole- 
doped materials should change sign for sufficiently high 
doping. 

Indeed, Ap was reported to be very different in 
electron-doped and hole-doped BaFe2As2 [27 , where it 
becomes vanishingly small for moderate hole-doping lev- 
els. Motivated by the prediction of Ref. [19], in this pa- 
per we report a systematic study of the in-plane resistiv- 
ity anisotropy for a broader doping range of underdoped 
Bai_a^Ka;Fe2As2 (BaK122), from x = (the parent com- 
pound) to compositions close to optimal doping, x = 0.34 
and Tc = 38 K. We find a change in the sign of the re- 
sistivity anisotropy for x > 0.235, confirming the fun- 
damental role of anisotropic spin-fluctuation scattering 
in the description of the resistivity anisotropy. Further- 
more, we also find that the temperature dependence of 
the anisotropy below T/v (i.e. in the magnetically ordered 
state) qualitatively changes for similar doping levels, with 
|Ap| monotonically increasing with cooling. This feature 
suggests that the magnetically reconstructed Fermi sur- 
face may undergo profound changes as positive charge 
carriers are introduced in the parent compound. 

EXPERIMENTAL 

Single crystals of BaK122 with sizes up to 7x7x1 mm^ 
were grown from FeAs flux as described in Ref. [28] . 
Potassium content in the samples was determined using 
electron probe microanalysis with wavelength dispersive 
spectroscopy (WDS), see [29 for details. Samples had 
typical dimensions of 0.5 mm wide, 2 to 3 mm long, and 
0.05 mm thick, and were cut from cleaved slabs along 
the tetragonal [110] direction (which becomes either the 
orthorhombic ao or bo axis below Tg). Optical imaging 
was performed on the samples while mounted on a cold 
finger in a continuous flow liquid helium cryostat (allow- 
ing precise temperature control in the 5K to 300K range) 
using a Leica DMLM microscope. 



Sample images in polarized light were used for visu- 
alization of the structural domains [30 and initial sam- 
ple selection was based on the clarity of domains in the 
image. Selected samples were mounted for four-probe 
electrical resistivity measurements, with contacts made 
by soldering 50 pm Ag wires using low-resistance Sn sol- 
dered contacts [31] . The first resistivity measurement on 
each sample was carried out using a flexible wire arrange- 
ment with no strain applied to the sample (free standing 
state). Samples were then mounted on a brass horse- 
shoe straining device, and strain was applied through 
the voltage contact wires by deformation of the horse- 
shoe, see Ref. [151 [2S] for the details of measurements. 
Strain in this configuration is applied along the tetrago- 
nal [110] axis, which then becomes the long orthorhom- 
bic ao axis upon cooling below T^. The current leads 
are arranged such that they do not apply additional 
strain. The strain was incrementally increased and, for 
each increment, temperature-dependent resistivity mea- 
surements were made and the domain structure imaged 
to determine the completeness of detwinning. When no 
structural domains were visible the sample was consid- 
ered detwinned. The strain needed to completely detwin 
the samples is estimated to be in the 10 to 20 MPa range 
based on our previous work [32]. Because of the small 
magnitude of the resistivity changes, detailed attention 
was given to optical imagaing to determine when the sam- 
ples were detwinned. 

RESULTS 
Polarized Microscopy 

In Fig. 1. we show polarized light images of the area 
between the potential contacts in BaK122 crystals used 
in this study, with x = 0.108 (non-superconducting), 
X = 0.202 (Tc = lOK), X = 0.235 (T^ = 26K), and 
X = 0.259 (Tc =28K). Images were taken on cooling at 
a temperature slightly above Tg and at base tempera- 
ture. The highest contrast is observed when the sample 
is aligned with the [110] tetragonal direction at 45° to the 
polarization direction of incident light (parallel and per- 
pendicular to the orthorhombic ao in different domains) . 
The contrast of domain images depends on the quality of 
the surface and the homogeneity of the samples. Domains 
are observed for samples with x ~ 0.26, but are not ob- 
served for samples with higher doping level, x = 0.34, 
which do not show any features in p{T) associated with 
the structural transition, see Fig. 2. 

Resistivity 

Temperature-dependent resistivity, see Fig. 2, was first 
measured in the samples without strain. In the twinned 
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FIG. 1. (Color Online) Polarized light microscopy of BaK122 
samples with doping levels as shown in the panels. Im- 
ages were taken at temperatures just above the tetragonal- 
to-orthorhombic structural transition (T > T^, right) and 5K 
(T < Ts, left), the latter showing formation of structural do- 
mains due to four different orientations of the ao and bo axes. 
The difference in color of the domains is due to spectral de- 
pendence of bireflectance, depending on the angle between 
the ao direction and the polarization plane of incident white 
light. Since bireflectance is proportional to the orthorhombic 
distortion, the contrast naturally vanishes in the tetragonal 
phase (right panels). 



crystal we may assume a roughly equal population of 
domains with the ao axis parallel to the current direction 
and domains with the bo axis parallel to it. Therefore, 
the resistivity in the free-standing crystals is close to the 
average value of the resistivities of the two orthorhombic 
axes, pt{T) = [pa(T) ^ pk(T)]/2. 

Our previous studies showed that when crystals are 
strained until no domains are visible in polarized mi- 
croscopy, about 90% of the bulk of the sample dis- 
plays the domain whose orthorhombic ao axis is oriented 
along the direction of the strain and, consequently, par- 
allel to the current. Therefore, in the detwinned state, 
we are predominantly measuring the ao axis resistivity. 
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FIG. 2. (Color Online) Temperature-dependent resistivity 
along the orthorhombic flQ-aXlS, Pai 

for samples of BaK122 
with five different potassium doping levels in the strain-free, 
twinned, state. The curves are shown using normalized plots, 
Pa (7") /pa (300ir ) . Arrows indicate the transition tempera- 
tures and were determined by the maximum in the derivative 
of the resistivity and from polarized optical microscopy. 



pa- Using this, we calculate the b axis resistivity from 
pb{T) = 2pt{T) — pa{T). The measured pa and cal- 
culated p5 are shown for samples with x = 0.202 and 
X = 0.235 in Fig. 3. Note that the application of strain 
has the opposite effect on the resistivity of the two sam- 
ples, with Pa < Pb for X = 0.202 (and all samples with 
smaller x), but pa > pb for samples with x = 0.235 and 
X — 0.259. In the bottom panel of Fig. 3, we plot the 
temperature-dependent anisotropy ratio pb / pa for all the 
samples studied, including those with x = 0. 

We summarize our results in Fig. 4. In the top panel, 
we combine the phase diagram determined from magneti- 
zation and neutron scattering measurements on polycrys- 
tals of BaK122 [33] and single crystals of BaCol22, with 
our measurements for the structural transition tempera- 
ture Ts and the superconducting transition temperature 
Tc. The former was determined from the maximum of 
the resistivity derivative dp/dT and from polarized mi- 
croscopy, while the latter was determined from the tran- 
sition mid-point. In the middle panel, we summarize 
the measurements of the in-plane resistivity anisotropy in 
electron- and hole-doped materials, using our current re- 
sults as well as previously published data. Specifically, we 
plot the maximum value of the temperature-dependent 
anisotropy for each case. It is clear that the data from 
different groups are in good semi-quantitative agreement, 
with minor differences coming most likely from poor con- 
trol of residual strain in the different experimental setups 
used for sample detwinning. The bottom panel displays 
the theoretical result of the model of Ref. [TF , and will 
be explained in the discussion section. 
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FIG. 3. (Color Online) Top panel. Normalized temperature- 
dependent resistivity, p(T)/p(300i^), for BaK122 samples 
X — 0.202 and x — 0.235. The red lines show the resistiv- 
ity along the ao-axis (pa), the blue lines show pb. The inset 
shows a zoom of the structural transition, with clear rever- 
sal of anisotropy from ph > pa for x = 0.202 to pb < pa for 
X = 0.235. Bottom panel shows the temperature-dependent 
anisotropy, pb / pa (T) , for samples with different K-doping lev- 
els. . 



Discussion 

The measured resistivity anisotropy of all electron- 
doped and parent compounds is such that Pb/Pa > 1 for 
all temperatures. In Bai_a;Ka^Fe2As2, this property per- 
sists for slight to moderate potassium hole-doping, rang- 
ing from p5/pa ^ 1-2 to 1.5, down to 1.02 for samples with 
X = 0.202. For samples with x = 0.235, this decreasing 
trend of pb / pa culminates in a sign-change of the resistiv- 
ity anisotropy Ap = p^^ — p^^ with ph/ pa ^ 0.956. Upon 
increasing doping even further, the anisotropy eventually 
vanishes due to the crossing into tetragonal phase. 

Since the sign-change happens above the magnetic 
transition temperature, as evidenced in Fig. 3, we first 
analyze our results without invoking the reconstruction 
of the Fermi surface due to long-range magnetic order. 
There are two main mechanisms that render the re- 
sistivity anisotropic in the paramagnetic/orthorhombic 
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FIG. 4. (Color Online) Top panel. Unified phase diagram 
of BaFe2As2 for the cases of electron (Co-substitution of Fe) 
and hole (K-substitution of Ba) dopings, showing the struc- 
tural, magnetic, and superconducting transition lines (after 
Refs. 33, 36 ). Data points (blue stars) are the measurements 
extracted in this study. Middle panel. Maximum in-plane 
anisotropy, pb/pa{T), as function of electron- and hole-doping. 
For both types of doping the anisotropy vanishes with sup- 
pression of the orthorhombic distortion close to optimal dop- 
ing, however, it reveals a steep maximum in electron-doped 
compounds (x=0.035 Co doping) versus a sign change in hole- 
doped materials with zero-crossing around x=0.17. Data 
points for the middle panel were taken from the following 
references: J. H. Chu et al. (black squares) Ref. 20 , M. 
A. Tanatar et al. (red circle) Ref. [15], and J. J. Ying et 
al. (green diamonds) Ref. [27 . Bottom panel. Theoreti- 
cal calculation of the maximum resistivity anisotropy Ap (in 
units of the residual resistivity due to impurity-scattering po) 
in the paramagnetic-nematic phase, as described in Ref. iTO] , 
as function of the chemical potential p (in units of the Fermi 
energy). Here, p > corresponds to electron-doping whereas 
p < corresponds to hole-doping. 



phase: orbital order (TTJ [12] and anisotropic scattering 
by spin- fluctuations [T9[ |34]. As pointed out by both 
Refs. [18] and [12], orbital order alone, as observed 
experimentally [16], gives the wrong sign of the resistivity 
anisotropy in the parent and electron-doped samples, in- 
dicating that the main contribution is coming from spin- 
fluctuation scattering. 
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In the model discussed in Ref. [19], the scattering 
by spin fluctuations becomes anisotropic below the or- 
thorhombic/nematic transition temperature because 
the fluctuations around one of the two possible stripe or- 
dering vectors, (tt, 0) and (0,7r), become stronger than 
the fluctuations around the other ordering vector. In 
this case, and when impurity scattering is also present, 
only electrons near the hot spots associated with the se- 
lected ordering vector are strongly scattered by spin fluc- 
tuations. Then, the sign of the resistivity anisotropy is 
determined by the projections along ao and ho of the 
Fermi velocity at the these hot spots. Roughly, if the 
hot-spot Fermi velocity has a larger component along 
the ho direction, electrons will be more strongly scat- 
tered by spin fluctuations when they travel along this 
direction, and > pa- We can then understand the 
sign-change in Ap = — as function of doping in a 
qualitative way: since doping changes the relative sizes of 
the hole and electron pockets of the Fermi surface, it also 
changes the position (and therefore the Fermi velocity) 
of the hot spots. In particular, electron-doping brings 
the hot spots closer to the ho axis, implying > pa^ 
while hole-doping brings the hot spots closer to the ao 
axis, yielding p5 < pa- To illustrate this general feature 
of the model, in the bottom panel of Fig. 4 we plot the 
calculated maximum anisotropy Ap as function of dop- 
ing in the nematic-paramagnetic phase, considering the 
particular case explored in Ref. [19 of critical magnetic 
fluctuations and large impurity scattering. 

A qualitative comparison with the experimental data 
in the middle panel of Fig. 4 shows that this model in- 
deed captures not only the sign-change of the resistiv- 
ity anisotropy, but also the asymmetry between electron- 
doping and hole-doping. However, the model does not 
seem to capture the peak of the maximum Ap at an 
intermediate electron-doping level. Since the theory is 
valid only in the paramagnetic state, for T > T/v, it sug- 
gests that this observed peak is related to the anisotropic 
Fermi-surface reconstruction in the magnetically-ordered 
state. Indeed, as pointed out in [14 , the doping level cor- 
responding to this peak in the maximum anisotropy coin- 
cides with a Lifshitz transition that was found in compar- 
ative ARPES/thermopower measurements [35 . In par- 
ticular, in this Lifshitz transition, a small hole-like recon- 
structed pocket disappears from the Fermi surface. 

This point sheds light on another interesting feature 
of our resistivity anisotropy measurements. As shown in 
the bottom panel of Fig. 3, the temperature-dependent 
anisotropy ratio Pb/Pa is maximum at for both the 
parent compound and the slightly hole-doped composi- 
tion X = 0.108. However, for the two other samples with 
higher doping levels, \pb — Pa\ increases monotonically as 
the temperature is decreased. One possibility, in view of 
the behavior observed in the electron-doped compounds, 
is that this change in the temperature-dependence of the 
resistivity anisotropy reflects a change in the topology 



of the reconstructed Fermi surface as function of hole- 
doping. Such a change should also impact other macro- 
scopic properties, and therefore could be verified by other 
experimental probes. On the theoretical side, it is inter- 
esting to note that Ref. [T8], investigating the Drude 
weights along ao and ho in the magnetically ordered 
phase, found different signs of the resistivity anisotropy 
depending on the properties of the reconstructed Fermi 
surface. 

Summarizing, we observe a change in the sign of the 
in-plane resistivity anisotropy when going from electron- 
doped to hole-doped BaFe2As2 compounds, in agreement 
with the general theoretical predictions of the anisotropic 
spin-fluctuation scattering model of Ref. [19 . The tem- 
perature dependence of the anisotropy below T/v also 
changes and becomes monotonically increasing with cool- 
ing, suggesting that changes in the reconstructed Fermi 
surface may also take place with hole doping. 
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